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2.0 INTRODUCTION 


The digital program TBAP has been developed to provide thermal 
boundaries in the D/M-relative velocity (D-V), dynamic pressure- 
relative velocity (q-V), and altitude-relative velocity (h-V) planes. 
These thermal boundaries are used to design and/or analyze Shuttle 
Orbiter entry trajectories. The TBAP has been used extensively in 
supporting the Flight Performance Branch of the NASA in evaluating 

candidate trajectories for the Thermal Protection System design 
trajectory. 
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3 *° discussion 

The method used in lew to establish t„ , 

« -Terence heat rate for a *• «» Tied 

reference heat rate i„t„ * . * J,t1tud e. input the 

«e definition of the ^ 

z:r : rrr ■-**' 

tlca] temperature input for that n* , 
end dynamic pressure are then found ‘hat panel. 0/M 

velocity point This “T-tpondinp to the altitude- 

s yuinc. rhis proces . , 

;r,::r • 

-~::r°7 ™ p are that ^ 

erature and that at™ = h ‘° dCterHi " e t«p- 

atmospheric density and free stream t 

can be modeled using the 1962 St „ ^erature 

9 ne 1962 Standard Atmosphere. 

The reference heat rate fir th 0 
<n Reference 1 is 

Vf = 0.62 (17700) v^-/ v \3.07 L 

W) (' -f rj ( 

where 

9 = stagnation point heat rate for a omp f * 

^ foot sphere (F\Tit/-r^ 2 / 

° ’ free *'-« density (s.ups/ftn 
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V * relative velocity magnitude (ft/sec), 

h w « surface enthalpy (BTU/1b m ), 

h w • - ?4 V t2) 

h Q * total enthalpy (BTU/lb m ), 

h = .24 T + V 2 (3) 

0 00 — 

50,063 

and T w = wall temperature (°R), 

T„ ■ 100o/Vf V 25 (4) 

M Wc / 

T = free stream temperature (°R), 

CO 


and e = surface emissivity. 

NOTE: Since T w , (equation (4)), is not independent of 0 ref > an initial 
value of T is assumed and tv/o passes are made through equation (1) 
in defining q fef . 


The individual panel or control point heat rate (Reference 1 and 2) 
for laminar, transitional and turbulent flow are defined by 


q , s n.„ , = A.C.q for P < R . . 

H panel M(lam) i vref oi 


hpanel *' °i(trans)’ q i(lo.ni) 
and 

hpanel ' °i(tu.'b) " q i(lam) c 



for P 0 . < R < K R q . 


5 . 


(5) 

( 6 ) 

(7) 


where A.. J deflection factor, 

C. = q i (la i) , for the i^ panel 
q ref 


: function of angle of attack, a , 
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R « Reynolds number behind a normal shock 


* transition onset Reynolds number for the 1th panel, 

* function of a, 

and y y « normal shock viscosity coefficient and a, b, c, and K are 

modeling constants. The deflection factor, A., (Equation (5)) equals 

one for surfaces not on the elevon or body flap and 

A. = 1.05 + (9.8 - a ) T 
1 5 

where a * angle of attack (degrees) 

T = .01965 - .0122, for 0 < 6 < 4 

T = .02766 - .0443, for 4 < 6 < 8 

T = .036765 - .11753, for 8 < 6 < 15 

T = .046696 - .2665, for 15 < 6 * 

and 6 * deflection angle (body flap or elevon) for surfaces on the 

body flap or elevon. 


C. and R . are stored in tables in TBAP as functions of a and interpolated 
1 01 

for the appropriate a. p is stored in a table. as a function of V and 
interpolated for the appropriate V. The surface temperature for the 
i ^ panel as defined in Reference 2 is given by 


T. - 1000 


panel \ .25 


.4761 c. 


v;h ere r. = surface enissivity for the 1^ panel. 

Using Equations 1 through 9, a rcgula falsi m-vthyd is used to solve for 
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the altitude at which panel temperature is equal to a critical 
temperature. 


That is 


•V, * N * [$] f (T c - T i> ( 

is initialized with an altitude h^, ^ calculated numerically by 

Ah . Vl Hid 

ST " T, - T 1-1 ’ 

and an iteration process continued until T, converges to T c within an 
input convergence criterion. 

The dynamic pressure is defined by 
q= \ PV 2 . 

In TRAP the free stream density, p, and temperature, T^ (Fouation 3) 
are computed using the SVDS (Reference 3) subroutine ATMOS which 
represents the 1962 Standard atmosphere. The drag acceleration is 

defined by 

D/M = ^ C D S . 


where = drag coefficient 
S = reference area 


and m - vehicl 


e mass. 


C-, is computed using the SVDS subroutine ARKC (Reference 3) 


i'f 
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4.0 INPUT DESCRIPTION 




The TBAP input data consists of data from the SVDS base data tape 
(Reference 3) for an entry run and data for the following variables 
which are in common arrays. 


VARIABLE 

TYPE * 

COMMON 

BLOCK 

COMMON 

LOCATION 

DESCRIPTION 

ALPHA 

R 

BLCK 

GWDAT1 (7S) 

Pitch angle of attack (RAD) 

MASS 

R 

BLCK 

GHDAT1 (128) 

Vehicle Mass (slugs) 

V 

R 

BLCK 

XLEC (507) 

Initial velocity (ft/sec) 

H 

R 

BLCK 

GNDATi (70) 

Geodetic altitude (ft) 

NC2M 

I 

GEN2 

W0RK2 (330) 

Convergence interval for 
temperature iterator 

NVAL 

I 

GEN2 

W0RK2 (331) 

Number of values for angle- 
of-attack table 

DT 

R 

GEN2 

W0RK2 (107) 

Specified tolarance factor 
for temperature iterator 

DV 

R 

GEN2 

W0RK2 (105) 

Increment velocity (ft/sec) 

NP 

I 

GEN2 

W0RK2 (108) 

TPS panel or control point 
number 

IOP 

I 

GEM? 

1/0RK2 (110) 

Option for atmospheric model 





= 0 1962 U.S. Standard 
atmosphere 





= 1 July atmosphere at 30 Deg 
North Latitude 





~ 2 Jan. atmosphere at 30 Deg 
North Latitude 





= 3 July atmosphere at 60 Deg 
North Latitude 





= 4 Jan. atmosphere at 6^ Deg 
North Latitude 

VM 

R 

GLN2 

\mY2 (104) 

Maximum velocity (ft/ sec) 






VARIABLE TYPE * 


COMMON 

BLOCK 


COMMON 

LOCATION 


IALP 

I 

GEN2 

W0RK2 (111) 

TC 

R 

GEN2 

W0RK2 (112) 

VALT (I), 
1=1, NVAL 

R 

GEN2 

W0RK2 (150) 

VALP (I), 
1=1, NVAL 

R 

GEN2 

V.'0RK2 (114) 

NBAP 

I 

GEN2 

W0RK2 (200) 

TSTOP 

I 

GNDAT2 

GMDAT2 (2) 

TMAX 

I 

GNDAT2 

GHDAT2 (3) 

I BAP 

I 

GEN2 

W0RK2’ (199) 


♦MOTE: Type "R" denotes REAL. 

Type "I" denotes INTEGER. 
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DESCRIPTION 


Option for angle-of-attack 
* 0 ramped schedule 
= 1 constant schedule 


Crltlal temperature (DEG) 
Velocity Table (ft/sec) 


Angle-of-attack Table (DEG) 


Flag determining thermal 
boundary option 

= 0 Does not call T8AP 

= 1 Calls TBAP 


Maximum Phase Time (sec) 
Maximum Case Time (sec) 


Flag for TBAP driver option 
= 0 Do not use TBAP driver 
* 1 Use TBAP driver 
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5.0 OUTPUT DESCRIPTION 


Figure 1 gives an example of the Intlallzatlon data which Is output 
from TBAP. Table I defines the data In Figure 1. 


Figure 2 gives an example of line printer output of thermal boundary 
data for a given panel from TBAP. Table II defines the data in 
Figure 2. 


r ' 


i 


Plots of thermal data from TBAP (see Appendix A Figures 3, 4, 5) are 
presented in the D/M-relative velocity, altitude-relative velocity, 
and dynamic pressure-relative velocity planes respectively for control 
point 2 (body flap). 





TABLE I. 

INITIALIZATION DATA 
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Input Identification 

Maximum number of iterations on critical temperature 
Number of values for angle-of -attack table 


Row 3-4 


Input identification 
REF AREA Reference Area (ft 2 ) 


WEIGHT 

VSTART 

VSTOP 


DELTA V 


Weight (lbs) 

Initial velocity (ft/sec) 

Final velocity (ft/sec) 

Velocity Increment (ft/sec) 

Specified tolarance factor for temperature iterator ( C F) 


Row 5-6 

Input identification 

CG (1) Center of mass X-component (in ) 

CG (2) Center of mass Y-coinponent (in) 

CG (3) Center of mass Z-component (in) 


■J 
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TABLE II 

THERMAL BOUNDARY DATA 

Row 1-2 

Input Identification 

NP TPS panel or control point number 

IOP Option for atmospheric model 

=0 1962 U.S. Standard atmosphere 
= 1 July atmosphere at 30 degrees north latitude 

= 2 Jan. atmosphere at 30 degrees north latitude 

= 3 July atmosphere at 60 degrees north latitude 

« 4 Jan. atmosphere at 60 degrees north latitude 

IALP Option for angle-of-attack 

5 = 0 ramped schedule 
= 1 constant schedule 
TC Critical temperature (°F) 


Rov/ 3-34 

Output identification 


REL VEL 

Relative velocity (ft/sec) 

ALTITUDE 

Altitude (ft) 

D/M 

Drag acceleration (it/sec ? ) 

OBAR 

Dynamic Pressure ( 1 Ms/ ft* ) 

DENSITY 

Density (slig/Yi 1 ) 

REY NO 

Not. ml shock Reynolds numk: 

MACH 

f'^ch number (NO) 

TEMP 

Surface temperature (°F) 


I 


i 



i 


i 


I 


i 1 

l 3 
« 

J ■ 
t ; 

i J 

1 ■: 

f | 


’* 






FLOW 


TABLE II (CONT'D) 
THERMAL BOUNDARY DATA 

Row 3 - 34 (cont'd) 

Output Identification^ 

Flow characteristic 
= 1 Laminar 

* 2 Transitional 

• 3 Turbulent 
Body flap setting (DEG) 

Elevon setting (DEG) 
Angle-of-attack (DEG) 


BDF 

ELV 

ALP 
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6.0 SUBROUTINE DOCUMENTATION 

The following sections define TBAP subroutines. The subroutines 
ATMOS and AR140C are current SVDS (Reference 3) routines and hence 
are not documented In this note. 

6.1 INTGI 2 

6.1.1 PURPOSE : INTGI2 

Subroutine INTGI2 calls the BAP routine. 

C.1.2 INPUT 

NBAP Flag determining thermal boundary option 
= 0 Does not call TBAP 
* 1 Calls TBAP 

6.1.3 OUTPUT 

None applicable to BAP 

6.1.4 ALGORITHM 

None applicable to BAP 

6.1.5 CALLING SEQUENCE 

Call INTGI2 
6.T .6 CONSTANT REQUIRED 
None required 

6.1. 7 SU BROUTINE REQUIRED 
BAP 

Other routines not applicable to BAP 
* r.'OTF: Only the portion of INTGI2 applicable to TRAP is documented 
in this design note. 


, J 




6.1.8 F10V1CHA R1 


srMT 




[cAU^P.t?. J 

3 L— 

[ r'.fJ ’-'LilJ 

ft'’ -* 

|^£/V7 3 MU { ! 


LCif -••*>{. 1- 

JV|16-1Z 

- ^r.a.^e„ 

I — * — 1 
^£7_7MN__1 


BAP, 


The remainder of subroutln 
INTGI2 is unchanged. 

(see Reference 3) 
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6.2 BAP 

6.2.1 PURPOSE : BAP Is the executive routine used to sequence other 

routines In calculating thermal boundaries. 

6.2.2 INPUT: 

NP TPS panel or control point number 

Y Initial velocity (ft/sec) 

VM Maximum velocity (ft/sec) 

DV Velocity Increment (ft/sec) 

I ALP Angle-of-attack option flag 

.* 0 ramped schedule 
= 1 constant schedule 

DT Temperature convergence criterion (°F) 

NC2M Cor 'ergence Interval specification 

6.2.3 OUTPUT : 

Input Identification labels 
QBAR Dynamic Pressure (lbs/ft 1 ) 

MACH Mach number (ND) 

H Altitude (ft) 

DOM Drag acceleration (ft/sec* ) 

V Relative velocity (ft/sec) 

Write LNDF1LE on UNIT 8 

6.2.4 ALGOR ITHM 

j 

See flowchart Sec. 6.2.8 j 

6.2.5 CA LLIN G SFOyENCF 

Call BAP 

. 

5 

i 

, 



6.Z.G 
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6.2.8 FLOWCHART (CONTINUED) 


DN No.: 1 .4-4-5 

Page: 21 








1 


I 


s 


4 


6.3 IPS 

6.3.1 PURPOSE: TPS 

Subroutine TPS calculates normal shock Reynolds 
number and the stagnation point convective heat 
rate for a one foot radius sphere. 


6.3.2 INPUT: 

IREN 

KNTTPS 

NPANEL 
RHOSL 
ANSVJ (3) 

IV 

MTPS 


NTPS 

VT 

T400K 

TEH? 

I i.M’ 


Reynold's number computational option flag 

Number of Integration step sizes since the 
last TPS Model Computational update 

TPS oanels and control points 


Density at sea level (slug/ft 5 ) 

Ratio for density at altitude to density at 
sea level 

Vehicle number 

Thermal Protection Option 


= 0 TPS is not simulated 

= 1 Initialization pass of TPS 

, 2 Execute TPS logic (Initialization completed) 

= 3 Print TPS summary 


Number of SVDS integration cycles per TPS 
integration step 


Relative Velocity (ft/sec) 


Time at altitude of 400000 feets 
Array of current atmosphere data 


Option for TRAP driver 

•s o Do not use 1BAP driver 
= 1 Use TPAP driver 


i 

J 
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6.3.3 OUTPUT ; 

QREF Stagnation point heat rate (BTU/ft*-sec) 
REN Reynolds number behind a normal shock 

6.3.4 ALGOR I THM 

See flowchart Sec. 6.3.8 

6.3.5 CALLING SEQUENCE 

Call TPS 


L 


PI 
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6.3.6 CONSTANTS REQUIRED 


Relative velocity table corresponding to VISC array 
DATA vt S/600* • lOOO* . 1500* ~200 q* *2500 • »3o00r,3b00. .4000. * 4 500. » 
t 5000* .5500* ,6?C0. ,650 q. ,7 0 00* , 7500 8C00 • 1 6500 . ,9000*, 9500. , 

l ICC0O«»tO500»|l 1 000 . t | 1 50« • * l 2000 • » 12500* • 1 30o0. *13500* » 14000. 

) , 1 MbOO. , 15GC0* , 15500. » i 6000* . t 6500 . . * 7000 * * l 7 5oO . . UOQO. * 

< 2pC0C. » 2 2CCC. ,25000. . ?550n. ,25 800 . , 2 6C0Q * * 260b 0 . , 26 1 00 • , 

> 3no00./ 

Coefficient of viscosity VS relative velocity 

DAT A V 1SC/.310E-6. .35IE-6 , .39<4E>6, .460E-6 , .54 IE-6 , .646E-6 , . 7 45E-6 . 

1 •850E-6 l *951E-6»*l065E-5».ll7i:-5».l28fc.-5»*13BE-5,.|47E-5, 

2 .|565E-r»».|655E-5».l7 3E>E-5,.l0lE-5,.ie65E~5,*l92t-5,*lV65E“5* 

3 .2L-6 ,.?0 3e-5, *2ICE-5 t .2J 7E-5 , . 22 7511-5 , . 2 40t-5 , .25 | E>5 , 

4 .26lt-5,*26Br-5,.275^.5».26?E-5».289E-5,*2925E-5i*297E-5, 

3 .3L-5,.3UE-5.*323E-5,.3 4 5E-5,.349E-5,.34bE-5,.34lE-5,.337E-5. 

6 .286E-5,*266f-6/ 


6.3.7 






6.3.8 FLOWCHART 
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6.4 HTRATE 

6.4.1 PURPOSE : HTRATE 

Subroutine HTRATE determines the flow characteristic 
(laminar, transitional, or turbulence) for a given 
panel and computes the corresponding heat rate. 

6.4.2 INPUT : 

ALPHA Angle-of -attack (RAD) 

NPANEL Numberof TPS panels and control points 
QDOTC Stagnation convective heat rate (BTU/ft 2 -sec) 

IV Vehicle number 

REY Reynolds number behind a normal shock 

DELVTR Angle of deflection (body flap and eleven) (RAD) 

IBAP Option for TBAP driver 

= 0 Do not use TBAP driver 
= 1 Use TBAP driver 

HP TPS panel or control point number 

6.4.3 OUTPUT : 

QDT Heat rate (BTU/ft 2 -sec) 

IPS Flov/ characteristic 

* 1 (laminar flow) 
c 2 (transitional flow) 

- 3 (turbulence flow) 

HENS Reynolds number et transition onset 

6.4.4 AUV.'RITM'* 

See flowchart Sec. 6.4.9 

6.4.5 CAt LJJjf. 

Call f'TRATF. 


’JkWfr 
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6.4 HTRATE 

6.4. 1 PURPOSE: HTRATE 

Subroutine HTRATE determines the flow characteristic 
(laminar, transitional, or turbulence) for a given 
panel and computes the corresponding heat rate. 

6.4.2 INPUT : 

ALPHA Angle-of -attack (RAD) 

NPANEL Number of TPS panels and control points 

QDOTC Stagnation convective heat rate (BTU/ft 2 -sec) 

IV Vehicle number 

REY Reynolds number behind a normal shock 

DELVTR Angle of deflection (body flap and elevon) (RAD) 

I BAP Option for TBAP driver 

s 0 Do not use TBAP driver 

* 1 Use TBAP driver 

NP TPS panel or control point number 

6.4.3 OUTPUT : 

QDT Heat rate (BTU/ft 2 -sec) 

IFS Flow characteristic 

= 1 (laminar flow) 

* 2 (transitional flow) 

- 3 (turbulence flow) 

REMS Reynolds number at transition onset 

6.4.4 ALG ORITHM 

See flowchart See. 6.4.9 

6.4.5 CALL I MG SE QUENCE 


Call HTRATE 


V RODUCIBILITY OP THE 
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Data (CU,K),K«I,|3) /.55,0, 


»89 0 0, 


data < c < 2 , k ) , k * i , t 3 > / • » s *4 o » 


23nC 


Data (C (3,k ) »k*i , 1 3 > /•c^hO* 


• • I 6 1 C » 

DATA (C <8 »K > »K»1 . |3) /.c2oO» 

• • I s |G» 


data (c<5,k» *k*i , i 3 » /. 1329 c. 


1 IqOi 


data ic <6,k i ,k=i , i3> /. 029 c. 


• P 7 80 1 


data < 0 ( 7 ,k i *k* i , , 3 ) /. 6 C 2 C 


23 O 0 


» • 


data (C(8,ki,k=i,i31 /. 2 M 35 , 


• 2 ^ 82 I 


Data <c<9,k).k=i,i3) /.o9 7 o. 


' *2280. 

data <cUo,K>,Kei,i3)/.oi 7 o» 
► • 1 1 n 0 « 
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• S I 70 . 
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•2b°9 , 
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•22CC, .2200. .2200 
.2993, .2516, ,2535 
•2617, .2637, .2660 
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, .2/00/ 
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Alpha table corresponding to LAMINAR QDOTL/QDOT STAG 

CUITA ALP/I0»»l 5 * *20.,25.» 27. 5,30* » 32.5,35., 37,5.**0«. ‘‘2.5. .50./ 

At transition onset VS Alpha and Panel number 

► I5COO./ , . tu i . fc 53nn, .HB2C0. .H08Q0. »3‘‘* t 00. . 
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. O0-,0..H 2 IC..3s‘>O..302n.. 2670.. 2^30.. 2300.. 23»0./ 
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( '< 1 R A 

• l 

3 * i 
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x ' A 

l • T ^ l 
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• 
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1 \ l.;! 1 
/ 1 t \ 
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■Cl*, LI ,L*I ♦U>/lM'r > rn03./ 

:S ‘ l V ,1. ) ,L* t , 1 ** )/3 9 70.»33r0. .2*50. ,2 3S0. , 
? 6?. , *950. ,8?>C» »7t;0« »720» ,A6o»/ 

■= ( ?a ,L ) |L= t , 1 *• ) / 1 4 > 1 J VCC». / 

S t r * ,H .LM , )H ,/i‘l. 5000*/ 

.S«22,i.l.L»l , l l » ) .MM. 13b: JO*/ 

. : 2 3 , L > . L = l • 1 *' ) > 1 *’ « 6 7 0 J W . / 

( 2“ ,1 ) , l.M , 1 ‘‘ ) / 1 ‘U 77' 1 ;iO >/ 

■i S 1 Jt S I L > I L r t I 1 ,f ) / *' l ' r . 0 • I l i ^ 7 0 t f 3 S C 0 • i 7 ? 7 0 • 

, »36?. %\?'' r • • : ICO. , IC:C'5. ,9/0. » tieo,/ 


. , » 7 5 c * . I - r» ; . , 9 i r . . r. • . , 7 HC . . 7 1 2 . , I'l "> • / 

• f ' { ■) ‘1 I 1 t L * 1 • 1 4 l / 1 4 » • 1 0 E '! 0 0 • / 

■- 2 V ,1 ) ,1 =1 , 1 '■ ) /t * 5870 . , 5 3 C 0 . • H63 °. ’ 

9 . 2 ^ 7 ?. ,;* 32 C. . 7173 . 1933 ./ 

iMf. co r r f " 7 nrfii ri'i t f : Reynolds nu-. l, i» at transition onset. 
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6.5 n 

6.5.1 PURPOSE: TI 


Subroutine TI computes the surface temperature as 
a function of heat rate. 


6.5.2 INPUT: 


QDT 

Heat rate (BTU/ft 2 -sec) 

NPANEL 

Array of TPS panels and control points 

NP 

TPS panel or control point number 

IV 

Vehicle number 


6.5.3 OUTPUT : 

Com; non 

TARY Array of panel surface temperatures (°F) 

6.5.4 ALGORITHM 

See flowchart Sec. 6.5.8 

6.5.5 CALLING SEQUENCE 

None 

6.5.6 CONSTANTS REQUIRED 

C = .24710577 E + 13 

6.5.7 SUBROUT INE REQUIR ED 

None 





■ V i iiiMiiiai^toi m m 
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m 


table 

6 - 6 -' purpose., table 

Subroutine TABLE u . » 

6-6 ' ? interpolate 

NVAL 

'<» ■ I . W al ,T r 0f va,ues 

m. , . , IMl „ ePen<fent *rra y 

vi flL De P en ^ent varfAhi 

XI va nable array 


XI 

6 * 6 -3 OUTPUT 

Jndepend) 

YJ 

6 * 6 * 4 ALGORITHM 

Dependent 


See f] owchart Sec * * 

6 - 6 - 5 ^ 

Ca " ™BLE (N Wl „ . 

6 ' 6 - 6 * x * •* vjj 

K = ) 

6 ‘ 6 ' 7 ^^WfiUIREo 

Hone 
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1.0 GENERATE INPUT PROGRAM (GIP) 

GIP is an auxiliary program which accepts a tape of thermal 
boundaries data from TBAP and plots the thermal boundaries using 
the TRWPLT (Reference 4) routine. The GIP was written to facilitate 
the method of inputting data into the TRWPLT routine. An example 
of the printed output from GIP is given in Figure 1. Examples of 
thermal boundaries which were plotted using GIP and TRWPLT are 
presented in Figures 2, 3, and 4. 

1.1 PURPOSE : 

GIP accepts a data tape from TBAP, processes the data, and outputs 
calcomp or microfilm plots. 


1.2 INPUT: 


iccomp 


KU.TIT 


NTRAII 


NCHmA 

Cr.'.jLT 
xla:,:i 
1l-". ‘ T 


Type plot indicator 

= 0 Calcomp 
= 1 Microfilm 

Data tape input unit 

Tape type selector 

- 0 Input data tape v/as generated using FORTRAN 
'..rife staler .amts • 

v. 1 input de t-?. tape was generated using NTRANS 
write statements. 

Plot sy.d.ol selector 

f !u It i pi i c a t i ve fa c tor 

Independent axis identification 

B •.‘•••oodon t axis ireutificulicn 




1 1 


* 




‘3 
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1.2 INPUT: (continued) 


i 


f 

i 

i 



r 




TITLE Graph title 

IPNF Total number of files on data tape for a given plot 

JANF Actual file number 

NVAR2 Dependent variable 

= 2 Altitude (ft/sec) 

= 3 Drag acceleration (ft/sec 2 ) 

= 4 Dynamic pressure (lbs/ft 2 ) 

L Number of cards to specify scale parameters 

M Humber of cards to specify axis identification 

1.3 OUTP UT: 

ICCOMP 
KUHIT 
NTRAN 
NCHAR 
CMULT 
XLABEL 
YLABEL 
TITLE 
I PRINT 


hOICEC 




Calcomp/microfilm indicator 
Data tape input unit 
Tape type selector 
Plot symbol selector 
Multiplicative factor 
Independent axis identification 
Dependent axis identification 
Plot title 
Print indicator 

= -1 Print plot titles only 

- 0 Suppress print 

- 1 Print plot titles and data points 

* n Print titles and first n data points (it ? 1) 

Data tape format indicator 

* 0 Tare contains r.rre Iran one type of record. 

- 1 Tape ecu '..jins cr.ly one type rf r f t< :rd. 




1.3 OUTPUT : (continued) 
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UNGER 


NOADV 


ENDLST 

ENDPLT 

ENDRUN 

REWIND 


End-of-file frame advance Indicator 

* 0 Frame or graph will be advanced normally. 

* 1 The variables specified by the first PLOT card 

following the ENDFIL card will be plotted on the 
current plot. 

Number of end-of-file marks to forward positioned over 
Frame advance indicator 

= 0 Frame or graph will be advanced normally. 

= 1 The variable in the PLOT list will be plotted 
on the same graph as the previous plots. 

Specifies variable to be plotted 

Plot list termination symbol 

Plot input termination symbol 

Job termination symbol 

Rewind the input data tape 


1.4 CALLING sequence 
N ot applicable 


1.5 constmts .required 
MVAR l =- 1 


l.o sit-routine REOuiRra 




1.7 FLOWCHART 
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FIGURE 1. Example of GIP Printer Output 
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4- ... i ; . : 

; i 1 i i i . i 
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ill' 

t t* -*-* * t 


i.. BODY FLAP 1 ; i ! ! 1 L . 1. I ..j . i • . ;. < I ; . ; 

i Critical Temperature * 2500°.! j j, * , j j j 1. _ . J. . 

I Deflection Angle « 16.3° i , | i I , | | ; _ . , • 1. 
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